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Abstract 


We present a numerical method for solving linear and nonlinear frac- 
tional partial differential equations (FPDEs) with variable coefficients. The 
main aim of the proposed method is to introduce an orthogonal basis of two- 
dimensional fractional Muntz—Legendre polynomials. By using these polyno- 
mials, we approximate the unknown functions. Furthermore, an operational 
matrix of fractional derivative in the Caputo sense is provided for computing 
the fractional derivatives. The proposed approximation together with the 
Tau method reduces the solution of the FPDEs to the solution of a system 
of algebraic equations. Finally, to show the validity and accuracy of the 
presented method, we give some numerical examples. 
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1 Introduction 


It is well known that fractional ordinary differential equations and fractional 
partial differential equations are important parts of the numerical analysis, 
statistical mechanics, and physics in recent years (see [5, 11]). 


Recently, several numerical techniques have been proposed by researchers 
for solving some classes of the fractional ordinary differential equations and 
the fractional partial differential equations (FPDEs). For example, Nagy [9] 
used the Sinc-Chebyshev collocation method for fractional nonlinear Klein— 
Gordon equations. Yin et al. [18] applied two-dimensional fractional-order 
Legendre functions for solving the numerical solution of the fractional partial 
differential equations. Chen, Sun, and Liu [3] used a generalized fractional- 
order Legendre function for solving FPDEs with variable coefficients. For 
more numerical methods, the interested reader may refer to [6, 4, 2, 8, 1, 13]. 


In the present paper, we apply two-dimensional fractional Muntz—Legendre 
polynomials (2D-FMLPs) to find a numerical solution of linear and nonlinear 
FPDEs with variable coefficients. 


This article is organized as follows. Review of Caputo fractional derivative 
is briefly provided in Section 2, and also we review one-dimensional fractional 
Muntz-—Legendre polynomials in this section. In Section 3, two-dimensional 
FMLPs are introduced. In Section 4, we introduce an operational matrix of 
fractional derivative in the Caputo sense. In Section 5, we approximate the 
solution of FPDEs, using the Tau method based on 2D-FMLPs. Section 6 
presents some examples of FPDEs to show the efficiency and the accuracy of 
the proposed method. Finally, a conclusion is expressed in Section 7. 


2 Preliminaries 


In this section, we review definitions and preliminary qualities of the frac- 
tional calculus, which are used in this paper. 


Definition 1. The Caputo fractional derivative of order a is defined as 


1 t u™ (7) 
Deu(t) = 4 Te) f Gonerdt, n-l<a<n, neéN, 
Pu) a=n,t>0. 


It can easily be shown that 
D°C =0, 


when C is a constant and 
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py = {° v € No and v < [al, (1) 
r 


rete, v € No and v > [a] or v ¢ No and v > [al], 


where [a] is the integer part of the positive real number a and No = 
£0,1,,..-}. 
Moreover, we have 


OD a;u,(t)) = ~ a;D°u;(t), 


where a;,i =1,...,n are constants. 


The fraction Muntz—Legendre polynomials L;(t;a@) on the interval (0, 1] 
are given by the following formula [4]: 


4 _4\i-k tel 
igo) =) Cal, Cy= “i [[(@t+ve+2). (2) 
k=0 . * v=0 


Remark 1. According to (2), the analytic form of L,(t;a@) can be written 
as follows: 


Loa) tie, (3) 
k=0 

where Cyr tk+a) 

aye ed 
bk i = x = T # (4) 

ki — k)IT(= +k) 

Also, we have 

Li(t;a) = PO =P (2#* 1), a > 0, (5) 
where plot ) are the Jacobi polynomial with parameters a, > —1; see 


[4, 16]. Using (5) and the recurrence relation between the Jacobi polynomials 
[16], we can obtain the following recurrence formula: 


Li+i (tsa) = af Li(t; a) — 6 Li-1(t; @), t=1,2 


ee) 


where 


ao - Cit al + 4 — VOQF 5 +12 - 1) - (- 17) 
‘ 2G4-1@ +2 )Gt+ o=1) 
fl i(¢@+2-1)Qi+24+) 
* @41)G4+ 4)2i1+ 4-1) 

1 


Lo(t;a)=1, Ln(t;a) = 2 +e 
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Figure 1: The graphs of L;(t; a) for various values of i and a 


Note that £;(1; a) = 1. Also, for a = 1, we obtain shift Legendre polynomials 
[15] as follows: 


L,(t;1) = P©°(2t-1), — t € [0,1]. 


For a 4 1, we have 
L,(t;a) = P&*—Y (24% — 1) ¢ PO (24% — 1), 
where PO) (o40 — 1) are fractional shift Legendre polynomials; see [6]. 


Remark 2. The FMLPs are orthogonal on the interval [0,1] with the or- 
thogonality relation: 


| tC a, (6) 
0 


~ Qi +1 
where 6,; is the Kronecker function. 
The graphs of FMLPs are shown in Figure 1. 


Definition 2. An arbitrary function u(t) that is integrable in [0,1], can 
be expanded as follows: 


u(t) = S°aLi(t;a), (7) 
i=0 


where ‘ 
a; = (2ia +1) | HOLLOW Gin 
0 


In practice, only the first (m+ 1)-terms of FMLPs are considered. Then 
it can be written by the infinite series of (7) as 
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m 


u(t) © um(t) = © aLi(t; a) = AT S(t; a), (8) 


i=0 
where 
= [ao,a1,... tial 3 
G(t;a) = [Lo(t;@), Li(tsa),...,Lm(tsa)]”. 
The following theorem shows that the approximation FMLPs converges to 


u(t). 


Theorem 1.{Error bound] Suppose D*“u(t) € C[0,1] for k = 0,1,...,m. 
If wtm(t) in (8) is the best approximation to u(t) from Mm = span{Lo(t; a), 
Ly(t;a),...,Dm(t;a)}, then 

Ma 
(ma +1)/2ma+1’ 


t) — m t w S 
llu(t) — Um (tlle SF 
where M, > |D™u(t)|, ¢ € [0, 1]. 


Proof. By applying generalized Taylor’s formula (see [3, 10]), we have 


tka ss 4: pmo a 
Also, 
m—1 pho M,pme 
t D*y(0t < 
Ia) » rast | =tgatD 


On the other hand, u(t) = A? (t; a) is the best approximation to u(t), and 


He ries D**u(0*) € Mm,o- So, it can be written as 


BR 


me tka ee 
u(t) ~um()I3 sll) = YY DM wlO*) 2 
k=0 
< Me i $2™ dt = Mz 
~T(ma + 1)? Jo T'(ma + 1)?(2ma + 1) 


Therefore, 


Ma 
(ma +1)/2ma+1 


llu(t) — tm (lle S 


Hence, the proof is complete. 
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3 Two-dimensional fractional Muntz—Legendre 


polynomials 


In this section, we define the fractional Muntz—Legendre polynomials in the 
domain 2 = [0,1] x [0,1]. To this end, we can define these polynomials as 
follows. 


Definition 3. Let {Z,(t;a)}92> be the one-dimensional fractional Muntz— 
Legendre polynomials on [0,1]. We call {Li(a;a)Lj;(t;8)}?5—9 the 2D- 
FMLPs (two-dimensional fractional Muntz—Legendre polynomials) on Q = 
(0, 1] x [0, 1]. 


Theorem 2. 2D-FMLPs are orthogonal on 2. 


Proof. For i # p and j # q, we have 
1 pl 
Jf tsle0) ty (ts 8) plese) Lal ts Baa 
0 Jo 
1 1 
=(f Lilasa)Lp(esa)de)( | £5(t:8)L4(t: Bat) =0, 
0 0 
and for 7 = p and j = q, we have 


i | [L;(a; a) ]2[L;(t; 8) Pdedt = ( | [Li(w, a) |dx)( i: (L(t; 8)[2at) 
ol 1 
~ Qiatl 2jB+1° 


Definition 4. An arbitrary function u(x,t) that is integrable in 2 = [0, 1] x 
[0, 1] can be expanded as follows: 


u(x,t) = D> So aij Li(a; a)L;(t; 6), (9) 


i=0 j=0 


where 
1 pl 
aij = (2ia+1)(2764+ yf | u(a, t) Li (x; a) L,(t; B)dadt. 
0 Jo 


In practice, only the first (m+ 1)(n+ 1)-terms of 2D-FMLPs are considered. 
Then it can be written by the infinite series of (9) as 


u(x,t) ¥ Umn(a,t) = 5° S > ass Li (a; 0) L;(t; B) = $7 (2; 0) ASCE; 8), (10) 


i=0 j=0 
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where 


(x; a) = [Lo(x; a), Li(2;@),...,Lm(a;a)]", 
$(t; B) = [Lo(t; 8), L1(t;B),..-, En (ts BY", 


and A is an (m+ 1) x (n +1) matrix as follows: 


ao0 401 *** Gon 
410 G11 °°" Gin 
Amd Ami *** Amn 


To calculate the coefficients a;; in (10), we use the Tau method, which is 
discussed in Section 6. 


Remark 3. The product matrix of FMLPs can be obtained as follows: 


(a; a)¢" (a; a)A ~ Ad(a; a), 


where A is an (m +1) x (m+ 1) product matrix for the matrix A. The 
elements of A are as follows, 


m 


Ajj = (23a +1) > angize, i,j7 =0,...,m, 
k=0 


where 


Gigjk =| L? (a) LF (x) Le (a) de. 


The following theorems state the existence of uniqueness and convergence 
analyses. 


Theorem 3. If the double series 


S- S- aig Li (x; a) L;(t; B) 
i=0 j=0 
converges uniformly to u(#,t) on the Q = [0,1] x [0,1], then 


ayy = (2ia+1)(276+ » f i u(a, t) Li (x; a) L,(t; B)dadt. 


Proof. Set 
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u(x,t) = S° > ais L;(x; a) L;(t; 8), 


i=0 j=0 


and let m,n be fixed. We can write 


=| | (2 Yi aejLalw;0)Lj(t; B)) Lm(a30) Ln (t; B)dedt 


i=0 j=0 
= Lat Liaia)En(aia)day( L;(t; B)Ln(t; B)dt) 
= a(f E3(e:a)de)( | L*(t; 8)dt) 
al 
= ai 5D) BH 


Lemma 1. If u(x,t) is a continuous function on [0, 1] x [0,1] and the double 
series 


S 2S. aij Li(a; a) L;(E; 8) (11) 
i=0 j=0 


converges uniformly to u(z,t), then the series (11) is the FMLPs expansion 
of u(x,t). 


Proof. The proof of this lemma is similar to that of [8, Lemma 1 ]. 


Lemma 2. If two continuous functions defined on [0,1] x [0,1] have the 
identical FMLPs expansions, then these two functions are identical. 


Proof. The lemma is proved similar to Lemma 2 in [8]. Oj 


Theorem 4. If the 2D-FMLPs expansion of a continuous function u(2,t) 
converges uniformly, then every 2D-FMLPs expansion converges to the func- 
tion u(x,t). 


Proof. It is a result of Lemmas 1 and 2. 


4 FMLPs operational matrix of fractional derivatives 


In this section, we introduce an operational matrix of fractional derivative. 
First, we present the following lemma. 
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Lemma 3. Let L;(t;a@) be an FMLP. Then the Caputo fractional derivative 
of L;(t;@) of order y > 0 can be obtained by the following relation: 


T(k 1 : 
= hang (ka + ) pka-y 
—¥+1) 


9 


where 


, . JO, ka € No and ka < 4, 
7 bei, ka € No and ka > ¥ or ka € No and ka > [9]. 


and by; is given in (4). 


Proof. By using relations (1) and (3), the proof of the lemma is clear. 


Now, the next theorem gives the operational matrix of fractional deriva- 
tive of FMLPs. 


Theorem 5. If d(x; @) is a fractional Muntz—Legendre vector given in (10), 
then 
D79(x;a) = D™ 9(2;), 


where D\ is the (m+1) x (n+1) operational matrix of fractional derivative 
of order y > 0. The elements of D are defined as follows: 


(y) host) a 
D3} (27 1) a : ; 


where 


» _ 40, ka € No and ka < 4, 
7 bki, ka € No and ka > ¥ or ka € No and ka > [9]. 


Proof. According to Lemma, 3, we have 


T(ka + 1) 


ka-y 
)= Dhar ka +1— au va 


On the other hand, the approximation of t**~7 by using (8) gives the follow- 
ing relation: 


pent ss S045. L5(t; a), (13) 


where 
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1 
Aj = (2ja+ yf tO-7 L(x; a)dt 
0 


I 1 
= (Qjat Nd, | prem ae de 
0 


s=0 
J 
1 
=(2ja+1 5 bs,; ————_—_—_-.. 14 
ie 2, “ka-y+satl1 a) 


Now, according to Equations (12)—(14), we can write 


T'(ka + 1) 1 
DI Li (ta) = (2ja+1) >) 5 ; bs. 50b 4 L;(t; a). 
par oee em T(ka+1—7) (ka—7y+sa+41) 
Thus, 
stat) 1 
D}, = (2ja +1) by 
oer > ahitthg—y41) (e+s)a-741 


k=0 s=0 


and this completes the proof. 


5 Approximate solution to the FPDEs 


In this section, we approximate the solution of the linear and nonlinear 
FPDEs by applying the operational matrix of fractional derivatives of 2D- 
FMLPs. 


5.1 Numerical solution of linear FPDEs with variable 
coefficients 


To study the linear FPDEs with variable coefficients, the following form is 
considered: 


O° u(x, t) O° u(a, t) _ 
a(t)—s a + Ot) aa C(x) +d(x)u(z,t) = f(z,t), (15) 


u(0, t) = go(t), u(z, 0) = ho(x), (16) 
where (a, t) € 2 = (0,1) x (0,1), 0 < a,6 < 1 and a, @ are the order of the 
fractional derivatives in the Caputo sense. In (15) and (16), the continuous 


functions a, b, c, d, go, ho, and f are known and the function u(z,t) is 
unknown. 


Ou(x,t) 
Ox 
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To solve problem (15)-(16), we approximate the known functions, the 


unknown function u(a,t), and the differentials —, — late by 
applying (10) and Theorem 5 as follows: 
u(a,t) =o" (#;0)Ud(t;8), f(x, t) = o" (a; a) F ot; 6), 17) 
a(z)~ AT g(a;a), —-(t) = b* (4; B)B, 18) 
cx) = C7 ¢(x;a), d(x) ~ E*¢(z;0), 19) 
ee 2 ~ 3 (a:a)(D)"U tt), 20) 
Mules 48 (e:a)(D)"VOlt: 8), 21) 
fal 
PED ~ 4 (a: aUD H(t; 8) 22) 


Substituting (17)—(22) in (15), we have 


$" (a; a) AT (D°)"UG(t; B) + $7 (a; (UD Bolt; 8) 
+ $7 (x; a)C7(D™)TUS(E; B) + $7 (@ a)ETUG(t; 8) = o7 (a; 0) FO(E; 8) 
or 
$7 (a; a)[A7(D*)7U + UD B + CF (D™)TU + ETU]¢(E; 8) (23) 
~ $7 (a; 0)F Olt; 8). (24) 


In addition, for the initial and boundary conditions (16), we have 


u(0,t) ~ o7 (0; 
u(x, 0) ~ 6" (a; 


US; B) ~ GGe(t 8), — go(t) ~ Go g(t; 8), (25) 


a)U9(0; 8) =o" (#;a)Ho, — ho(x) =o" (a;a)Ho. (26) 


Now, by using the Tau method in [14], we could generate m x n equations 
from (24) as follows: 
a at 
7 / L&(x)$" (a; a)[AT(D*)?U +UD B+ CT (DY)TU + ETUI A(t; A)LE (t)daxdt 
0 0 
‘ll 1 
~ ff tx@oTaayroua\Ez@ded. (7 
0 0 


By using the orthogonality property (6), we have 
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[ 0 


a L é 
[ff eee? esanotes 6948 deat = : | 
0 0 0 


1 
(iat) (@jB+1) 
0 


0 
0 


fori =0,1,...,m—1 andj =0,1,...,n—1. Therefore, we rewrite (27) as 
follows: 


A’ (pF 40D B.C (Dey ULE yak (28) 


Also, we generate 2n + m+ 1 equations from conditions (25) and (26) as 
follows: 


¢' (0;a)U~Go, j=0,1,...,n, (29) 
U¢(0; 8) ~ Ho, i=0,1,...,m—2. (30) 
Now, we have a system of the mn+m-+n-+1 algebraic equations of Equations 
(28)—(30) with the mn +m+n-+1 elements of the unknown matrix U. After 


solving this matrix system, by using the fixed point method, we can find the 
approximate solution Um.n(a,t) = 67 (x;a)U4(t; 6). 


5.2 Numerical solution of nonlinear FPDEs 


To study the various types of nonlinear FPDEs, the following forms are con- 
sidered: 


(a) Time-fractional convection-diffusion equation: 
A class of time-fractional convection-diffusion equation may be written as 
follows: 


Bye = tae — Cle + N(u) + f(a,t), (31) 


u(0, t) = go(t), u(1,t) = gi(t), u(z,0) =ho(x), O<ax<1,t>0, 
(32) 
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where N(u) is a nonlinear operator. In this paper, we choose one or more 
terms of the nonlinear operator N(u) in the form of u*,k CN, wuz, or UUge- 


(b) Time- fractional-order Klein—Gordon equation: 
A class of time-fractional order Klein—Gordon equation reads as follows: 


O° u(x, t) at 07 u(a, t) 
Ot Ox? 
u(0, t) = go(t), u(1, t) = m(t), u(ax, 0) _ ho(a). 


+ au(a,t) + bu? (x,t) + cu(a, t) = f(a,t), 


(c) Time- fractional-order gas dynamics equation: 
A class of time-fractional order gas dynamics equation reads as follows: 


O° u(a, t) Ou(x, t) _ 
aco a u(x, t)— —u(z,t)(1 — u(a,t)) = f(a,t), 
u(0,t) =go(t), — u(w, 0) = ho(z). 


To solve problem (a), (b), or (c), we approximate the unknown function 
u(x,t) as follows: 


u(x,t) & Um n(a,t) = S- So ug Lila) L;(t; a) 


i=0 j=0 
= UT Y (x,t) =U" (d(2) ® o(¢4)), 
where ® denotes the Kronecker product and L;(a”) = L;(x; 1) are the Muntz- 


Legendre polynomials with a = 1. Also ¢(x), ¢(t;a), U, and Y(z,t) read 
as 


o(2) = o(a; 1) = [Lo(x), Li(a), ea pe De Is 
o(t; a) = [Lo(t;@), Li(t;a),...,Ln(tya)]”, 


U00 Lo (a) Lo(t; a) 
uo1 Lo(«)L1(t; a) 
Uon Lo(x)Ln, ; a) 
w10 Li (x)Lo(t; a) 
u=| : |, Y= 
Ugj—1 Lj (x) L5—1(t; a) 
Ui,j Lj (x) L; (te) 
UWij+1 L(x) L541 (t; a) 
Um,n—-1 Lm(x)Ln—1(t; a) 
Umn L Lim («) Ln (t; a) 4 
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Furthermore, the first and second order derivative operator matrices of the 
vector T(#,t) can be obtained by 


a 

5g tet) = DeT(a,t), De = DO) Teer 
62 
aga tat) ~ D2Y(z,t), D2? =D? @ Ina, 


where D) and D®) are the first and second order derivative matrices of the 
vector ¢(), respectively. 

The Caputo fractional derivatives operator of order a > 0 of the vector 
Y(a,t) can be obtained by 


DPT («,t)~ DET (e,t), DS = Ins @ DY, 


where Dio is the Caputo fractional derivatives matrix of vector $(t; a), which 
is obtained in Theorem 5. 


Remark 4. The product matrix of FMLPs vectors can be obtained as 
follows: - 
Y(a,t)T7 (a,t)U ~ UY (a, t), 


where U is an (m+1)(n+1) x (m+1)(n+1) product matrix for the vector 
U as follows: 


f= [59 G2, Gm FLO Fm m0 Emmy 
where U?*4,p = 0,1,...,m,q=0,1,...,n, isan 1 x (m+1)(n+1) matrix as 
Gs — (084,088... ORS, 08,..., 084... 084, OFM 
and each element of U?4 is obtained as 
m n 
ORT = (2k+1)(2la+1) s S- UsGinekis Pek =O, 1,...59%, g,b= 0,100.40, 
i=0 j=0 


where gijpqki is given by 
a 
siinant = f | L,(a)L;(t; a) Lp(a) L(t; a) L(x) Li(t; a)dedt. 
0 Jo 


Remark 5. Assume that f (x,t) ~ F?Y(a,t) and g(a,t) ~ G7 Y(z,t), where 
F and G are an (m+ 1)(n + 1)-vector as 


T 
E SWp0s Fos Pes ond iios “a tie te legedeg tis 6 lament 3 


Tr 
G =[90,0, 90,1; +665 GJO,n> 91,09 ++ +5 Gi,j-1) Gi,j, Gi,j+1s-+- Grn tran ) 
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where, from Definition 4, each element of F' and G is obtained as 


fig =Qi+ NQja 1) ff He ntsaybs(tsa)dedt 


i? “al 
Gi,j =(21+ 1) (27a yf | g(a, t)L;(x)L,;(t; a)dadt. 
o Jo 
We can write 
f(a, t)g(2,t) ~ APY (a, t)Y7 (2, t)B ~ ATBY (x,t) ~ BT AY (a, t). 


Lemma 4. The nonlinear operator N(u(z,t)) = u*(x,t), where u(z,t) ~ 
UT Y (x,t), can be approximated by 


u® (x,t) ~ UTU®11Y(a, t). 
Proof. From Remark 5, we have 
u?(a,t) = u(a,t)u(x,t) ~ UTUY (2,2), 
u®(a,t) = u?(a, t)u*-? (a2, t) ~ UTUOY (a, t)u(a, t)u*-3(c, t) 
~ UTUY (a, t)Y? (2, t)Uu*-3(z, t) 
~ UTE? (a, thu*3(a,t) ~~ UTO 1 (2, 2). 


Remark 6. The nonlinear operator Ni (u, uz) = wuz and No(u, Ure) = UUae, 
can be approximated by 


Ni(u, Ug) = wu, ~ UTD,UY (z,t), 
No(U, Uex) = Wee ~ UF D2UY (z, t). 

Now, to solve one of the problems (a), (b), or (c), for instance, problem 
(a), with rhe nonlinear operator N(u) = wuz + UUz, + u?, by applying the 
above approximations, we rewrite (31) as 

[U7 (D? — D? + cD) — N]Y(z,t) = F7Y (2,2), (33) 
where 


N(u) & NY(z,t) =UT[D,U+D?0+U0)Y(z,t), — f(z,t) & FTY(z, 2). 


As in the Tau method in [14], we should generate (m — 1)n equations from 
(33) as follows: 


[U7 (D¢ — D? + cD) — N] | | Y (a, t)Li(x)L;(t; a)dardt (34) 
0 0 
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=" f x (x) L; (t; a)dax 
=F | i T(a, t) L(x) L;(t; a)dedt, 


fori = 0,1,...,m—2 and j = 0,1,...,2—1. On the other hand, by using 
the orthogonality property (6), for 7 = 0,1,...,m—2 andj =0,1,...,n—1, 
we have 


Oo 


oo 


[ : Y(a, t)Li(w) Lj (t; a) dxdt = 


(Qi+1)(2ja+1) 


0 
1 
0 


= es mere 


Therefore, we rewrite (m — 1)n equations (34) as follows: 
[U7 (D? — D2 +cD,) — N] = FT. (35) 
Also, we generate 2n + m+ 1 equations from conditions (32) as follows: 


O'@eOQ@hal=G, w@)=Chete), 7=O,1..4n, 
U* (6(1) ® In41) = GT, git = GT ¢(t; a), J = ,1,.-.,7, 
U" (Im+1 ® $(0;a)) = HG, ho(x) = Hp 4(2), i=0,1,...,m—2. 


<5 


Now, we have a system of the mn +m-+n-+1 nonlinear algebraic equations 
of (35)-(38) with the mn +m-+n-+4 1 elements of the unknown matrix U. 
After solving this linear system, by using the fixed point method, we can find 
the approximate solution Um n(z,t) =U? Y(z,t). 


6 Numerical illustration 


In this section, we present some examples of linear and nonlinear of FPDEs 
to show the efficiency of the proposed method. The results will be compared 
with the exact solutions and other methods. The accuracy of the presented 
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method is estimated by the absolute error Ey,» and the maximum absolute 
€ITOY €m,,, Which are given as follows: 


Em,n = max{|u(O;, Bj) — Um (4, By)I}, 


where 6; and (6; are defined as, 


(9:,8;) = (2,t8), t=0,1,...,m, j=0,1,...,2, 


where x; and t; are Chebyshev-Gauss—Lobatto points with the following 
relations: 


1 1 Tt Eig 
r= = — ~=cos—, ;= 7s —. 
2 2 m J 2 2 n 


Example 1. Consider the linear FPDE with variable coefficients 


d°u(a,t)  APu(a,t) _ 1 = 1 2 
te + ao =Fara* ‘tra-p™ © co 
u(0,t) = u(a,0) = 0, (a, t) € (0,1] x (0,1), O<a,6B <1. (40) 


The exact solution is u(x,t) = at. Applying the presented method to 


solve (39) and (40) with m =n = 2, for a = $ and § = 3, we get 


2.5000e — 01 2.0000e — 01 5.0000e — 02 4 ; oT 
U = | 2.0000e — 01 1.6000e — 01 4.0000e — 02} ~ 5 25 5K 
5.0000e — 02 4.0000e — 02 1.0000e — 02 30 35 66 
Therefore, the approximation solution is 
tr, 1 1 
ug,2(x,t) = b (2; Uo 5) 
a eo 1 
13 2 10x — 12 g |i 2 2 3Vt — 2 = 
30 35 Too) [10t—-12Vt+3 


which is actually the exact solution. The results are shown in Tables | and 2. 
Table 1 shows the comparison of the presented method with other methods 
in [3, 14]. From Table 1, it is seen that we obtained a lesser error. Table 2 
shows the maximum absolute error €,,,,, for various values of m,n,a, and £. 
It can be seen that the approximate solution is obtained from the proposed 
method is in a good agreement with the exact solution. 


Example 2. Consider the linear FPDE with variable coefficients 
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TEE Absolute error fr 


‘Absolute error for m=n=6 


Error 
Error 


TEE Absolute error for m=n=7| 


Error 


ie? TERE Absolute error for m=n=9 [MN Absolue err for m=n=10] 


Error 
Error 


Figure 2: The absolute error function with a = 8 = # and m = n = 5,6,7,8,9,10 for 
Example 2 


3 


4 


4S 


LOS cfm 


5 


5 10 15 20 5 10 15 20 


Figure 3: The logarithmic graphs of €m,m for Example 2 
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Table 1: The absolute error for example 1 
the presented method method in [3] method in [14] 
(x, t) 
vps I=2> 5 a= f=4 
m= n= m= ih = a= B = 

(0.1, 0.2) 3.4694e — 18 3.47e — 17 1.20e — 04 

(0.5, 0.2) 0 0 7.61e — 03 

(1.1, 0.2) 2.7756¢ — 17 5.55e — 17 2.11e — 02 

(1.5, 0.2) 5.551le — 17 1.39e — 16 3.27e — 02 

(1.9, 0.2) 0 1.67e — 16 4.65¢ — 02 

Table 2: Maximum absolute error for Example 1 
€55 4.3696e—-14 3.8858e—15 3.1086e — 15 
€7,7 7.1054e—15 1.4433e-—15 1.1102e — 15 
.  O%u(a, t) OP u(a, t) 
sintg—s a + (1 + cost) ae /cu(z,t) = g(2,t), 
u(0,t) = u(x, 0) = 0, (x,t) € (0,1] x (0,1), 0< a,6 <1, 
where 
l-a Qe2-% 942-8 3 

D=t si +(1 14 al 

g(a, t) sint(Fa— a Tara)’ (1-2) (x( aaa te ) 


The exact solution is u(x,t) = t?(a—«?). The results are shown in Table 
3 and Figures 2 and 3. Figure 2 shows the absolute error function Em,» for 


a=P= 


and various values of m and n. Figure 3 shows the logarithmic 


graphs of Em n (lOZ19 m,n) for various values of a. It is seen that the absolute 
error converges to zero by increasing m and n (m= 7). 


Table 3: Maximum absolute error for Example 2 


inn @=B8=1  a=S=" a=f=! a=p=5 

€5,5 1.6256e —03 7.4183e—04 8.1055e—05 3.6924e — 04 
£66 3.2465e —04 8.1570e—-—05 4.6443e—05 2.0615e — 04 
E7,7 4.4705e-—05 4.4833e-—06 2.6524e—05 1.3334e — 04 
Egg 4.3807e —05 2.1630e—14 1.6770e—05 9.0206e — 05 
£99 2.1148e —O7 1.1367e—13 1.0939e—05 6.3418e — 05 
€10,10 1.0694e—14 1.5205e—13 7.3022e—06 4.6034e — 05 


Example 3. Consider the following time-fractional order Klein—Gordon 


equation 
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Ti Exact solution for a=1/2, [Tia Approximate solution for m=n=4 & a=1/2, [BREE Absoitte error for mana Bat 


Figure 4: Comparison between the exact solution and the approximate solution with 
absolute error function E4,4(x,t) for a = 1/2 for Example 3 


O°%u(z,t)  6?u(z,t) - 2 * ss 
ote Ox? eh) I'(3 — a) a 
u(0,t) = t?, u(x, 0) = 2”, u(1,t) =1+#°, (a, t) € (0,1) x (0,1], 0<a<1. 


The exact solution is u(a,t) = 2? +¢#?. Applying the presented method 
to solve this problem for a = 4 with m = 4,n = 4, we get 


6.6667e — 01 3.8095e—O01 2.1429e—O01 6.3492e—02 7.9365e — 03 
3.8095e — 01 1.3292e —13 —9.7353e—14 3.7575e—14 5.2979e — 15 
U = |2.1429e —01 1.0596e —14 —7.848le — 15 —2.3027e — 14 —3.2506e — 15 
6.3492e — 02 —3.3587e — 14 2.689le— 14 —1.3973e — 14 —1.846le — 15 
7.9365e — 03 —1.025le— 14 8.1934e—15 1.3784e—15 5.2590e — 17 


Therefore, the approximation solution is obtained as 


ua a(t) = 88 (a; SUL 5) 


r 1 - 1 
1 1 
3a2 —2 3t2 —2 
1 1 
= 102 — 1272 +3 x U x 10¢ — 1242 +3 
all 3 ul 3 
30a 2 + 3522 — 60a —4 30¢2 4+ 35¢2 — 60t — 4 
3 1 3 1 
[2102 — 28022 — 60x2 +1262? +5 210t — 280¢2 — 60t2 + 126t? +5 
1 1 iL, 3 
= 0.000045a — 0.00006a2 — (3t2 — 2)((1.2094e — 13)x2 — (1.6947e — 12)x2 
+ (3.1530e — 14)a@ + (1.2916e — 12)x? — 3.8095e — 01) + #? + ((8.9304e — 14) 


3 A 
(7.9339e — 14)a2 — (3.6375 — 15)? + (6.6263e — 15)x” + 7.9365e — 03)(210t 
3 1. 1: 3 
— 280¢2 — 60t2 + 126t? + 5) + (30¢2 + 352 — 60t — 4)((8.9757e — 13)x 


2 


iw 


3 
(8.7501le — 13)a2 — (1.1285¢ — 13) 


(1.7368e — 13)2 
3 
2 


6.3492e — 02) 


+ (10t — 1242 + 3)((1.1724e — 13)x2 — (1.3530e — 12)x 
+ (2.8673 — 14)a + (1.0324e — 12)a? + 2.1429e — 01) + 0.33335. 
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TE Exact solution for a=173 
© Approximate solution for m=n=10 


x 
= 


Exact solution for a=2/3 
Approximate solution for 


Error 


Figure 5: Comparison between the exact solution and the approximate solution 


absolute error function Emn form =n =10 and a= 


GEE Absolute error for man=10 & oat 


TEE Absolute error fo &ant2 


TERE Absolute error for men=10 & a=2/3 


with 
5,0 = 4,0 = 2 for Example 3 
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Figure 4 shows the comparison between the exact solution and the approxi- 
mate solution with absolute error function E,,,,(z,t) = |u(z,t) — Un n(2, t)| 
for m = 4,n = 4, and a = 1/2. Moreover, Figure 5 shows the comparison 
between the exact solution and the approximate solution with absolute error 
function Ey for m = 10,n = 10 and a = 1/3,a = 1/2,a = 2/3. Also, It 
can be seen that there is a high level of accuracy. 


Example 4. Consider the nonlinear FPDE 


u(x,t) 
ot 
u(0, t) = Balt"), 


+ u? (a, t) + u(x,t) = g(z,t), 
u(z,0) =e”, (a,t) € (0,1] x (0, 1], 


0<a<l, 


where 
g(x, t) = e**(E,(-t*))’. 


The exact solution is u(a,t) = e” Ey(—t®), where E,(—t®) is the Mittag— 
Leffler function defined in [12]. 
The numerical results are shown in Figures 6 and 7. Figure 6 shows the 
comparison between the exact solution and the approximate solution with 
absolute error function Ey, for various values of m,n, and a. Also, Figure 
7 shows the logarithmic graphs of €m,» (log;9 €m,n) for various values of a. 
It is seen that the absolute error converges to zero by increasing m and n. 


Table 4: Maximum absolute error for various values of m,n and a for Example 5 


Em,n a=0.2 a=0.7 a=0.8 

£33 1.3809e +00 2.5202e —02 1.5926e — 02 
EAA 6.7523e —O01 2.2084e—03  1.1443e — 03 
£55 1.5944e —01 1.8012e—04 7.4059e — 05 
£6,6 5.5797e — 02 1.3632e—05 4.3405e — 06 
€7,7 1.6813e —02 9.6897e —O07 2.3238e — 07 
E88 4.9588e —03 6.925le—08 1.1452e — 08 
E99 1.5654e —03 4.4995e—09 5.2297e — 10 
E10,10 4.9928e—04 2.7060e-—10 2.2264e — 11 
E1111 1.8696e —04 1.5228e—11 8.9617e — 13 
E1212 4.0219e-—05 8.0180e—13  1.5720e — 13 
€13,13 1.1839e-—05 1.5720e—13  1.7497e — 13 
Ey414 3.2635e-—06 1.1945e—-13 4.4941e — 13 
€15,15 9.0875€—-O7 2.3048e—12 8.7752e —13 


Example 5. Consider the nonlinear time-fractional gas dynamics equation 


[17] 
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[EE Exact solution for a=172 4a 
@ Approximate solution for m=n=9 & a=1/2 % [EMAbsotute error for man: 


:xact solution for a=3/4 HERE Absolute error for m=n=12 & a=3/4 
‘Approximate solution for 


TERE Exact solution for a=t 
© Approximate solution for man=18 & a=t 


Figure 6: Comparison between the exact solution and the approximate solution with 
absolute error function Em,n for various values of m,n and a for Example 4 
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oS mm 


108s cF nm 
Log, 
LOG Fm 


O° u(x, t) 
ote 
u(0,t) = Eg(t®), u(a,0) =e", (x,t) € (0, 1] x (0, 1). 


The exact solution is u(x,t) = e~* E(t). 

Figure 8 shows the approximate solutions obtained by the proposed method 
for t = 1 and various values of a. Similar results were obtained by Tamsir 
and Srivastava [17], by employing the fractional reduced differential trans- 
form method. For a = 1, the exact solution is u(x,t) = e’~*; see [7]. Figure 
9 shows the comparison between the exact solution and the approximate so- 
lution with the absolute error function for a = 1. We obtain a high level 
of accuracy, similar results in [17, 7], and we see that these results confirm 
the theoretical results. Also Figures 10, 11, and 12 show the comparison 
between the exact solution and the approximate solution with the absolute 
error function for various values of m,n, and a = 0.2, 0.7,0.8. In addition, Ta- 
ble 4 shows the maximum absolute error €m,,, for various values of m,n, and 
a = 0.2, 0.7,0.8. Figure 13 shows the logarithmic graphs of €m » (Logio€m,n): 
It is seen that the absolute error converges to zero by increasing m and n. 


Example 6. Consider the inhomogeneous fractional convection-diffusion 
equation 


O%u(x,t) 0? u(z,t) a Ou(x,t) 
ote Ox? Ox 


+ u?(a,t) + u(z,t) = g(2,t), 0<a<l, 


where . 
g(x,t) =TO +a) + (2? +2°) 4+4°% +07 + Qe — 2. 


The exact solution is u(x,t) = 27 +t. 
The numerical results are shown in Figure 14. 
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Figure 8: The approximate solutions obtained by proposed method for t = 1 and various 


values of a for Example 5 


TER Exact solution for a=t 
© Approximate solution for m=n=10 & a=1 


UCR) Uy gol 8) 


Figure 9: Comparison between the exact solution and the approximate solution 


absolute error function F19,19 for a = 1 for Example 5 


TERE Absolute error for man=10 & a=1 


with 
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‘act solution fo 6 [i Absolute errr for m=n=10 & a=02 
Approximate solution for m. 


Error 


‘act solution for a=0.2 Absolute error for m=n=20 & a=0.2 
pproximate solution for m: 


saa) 


Error 


Figure 10: Comparison between the exact solution and the approximate solution with 
absolute error function Emn for m =n = 10,15, 20, and a = 0.2 for Example 5 
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Exact solution for a=0.7 
4 L el [Absolute enor for mana10 & w=0.7| 
‘© Approximate solution for m=n=10 x10" a 


x10" 
Exact solution fo TEE Absolute error for men=16 @ a=0.7 
Approximate solution fo 25, 


Error 


Figure 11: Comparison between the exact solution and the approximate solution with 
absolute error function Em.n form =n = 10,15, and a = 0.7 for Example 5 
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Absolute error for m: 


Exact solution for 0=0.8 
Approximate solution for 


Error 


B0=08 


xact solution for a=0.8 
pproximate solution for ad [ERR Absolute error for m: 


Error 


Figure 12: Comparison between the exact solution and the approximate solution with 
absolute error function Em,n form =n = 10,15, and a = 0.8 for Example 5 
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Figure 13: The logarithmic graphs of m,n for « = 0.2,0.7 and a = 0.8 for Example 5 
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TEE Absolute error for m=n=10& a=1/8 


[TE Bxcet solution for a=tva 
Approximate solution for m=n=10 & a=1/3 


TE Bxcet solution for a=1/2 
Approximate solution for m=n=10 & a=1/2 MER Absolute error for mon=10 8 a=1/2 


[Eye xact solution for a=a/a 
‘@ Approximate solution for m=n=10 & a=3/4 Tp bsolte ertor for 


10 & a=3/4 


Figure 14: Comparison between the exact solution and the approximate solution with 


absolute error function F10,190 at a = i, s and 3 for Example 6 
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7 Conclusion 


In this paper, we applied a basis of 2D-FMLPs to obtain the numerical so- 
lution of linear and nonlinear FPDEs with variable coefficients. To get the 
unknown coefficients 2D-FMLPs, we used the operational matrix of frac- 
tional derivatives of FMLPs together with the Tau method. The results of 
the numerical examples and the comparison with other methods showed the 
efficiency and accuracy of the proposed method. 
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